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Relativistic viscous hydrodynamic fits to RHIC data on the centrality dependence of multiplicity, 
transverse and elliptic flow for ^/s = 200 GeV Au+Au collisions are presented. For standard 
(Glauber-type) initial conditions, while data on the integrated elliptic flow coefficient V2 is consistent 
with a ratio of viscosity over entropy density up to rj/s ~ 0.16, data on minimum bias V2 seems to 
favor a much smaller viscosity over entropy ratio, below the bound from the AdS/CFT conjecture. 
Some caveats on this result are discussed. 



The success of ideal hydrodynamics for the description 
of heavy-ion collisions at the Relativistic Heavy- Ion Col- 
lider (RHIC) has led to the idea of a quark-gluon plasma 
behaving as a "perfect liquid" , with a very small ratio 
of viscosity over entropy density 0, [3, H, An an- 
swer to the question "How perfect is the fluid observed 
at RHIC?" can, however, not be found using ideal hy- 
drodynamics, but must involve a controlled quantitative 
understanding of non-idealities, e.g. viscous effects. If 
hydrodynamics can be applied to RHIC physics, then 
relativistic viscous hydrodynamics should be able to pro- 
vide such an understanding. In particular, if one has 
control over the initial conditions, it should be possible 
to determine the size of various hydrodynamic transport 
coefficients, such as the shear viscosity, by a best fit of 
viscous hydrodynamics (VH) to experimental data. In 
this Letter, we aim to take a step in this direction. 

For RHIC physics, since particle number in the quark- 
gluon plasma is ill-defined, the relevant dimensionless pa- 
rameter for VH is the ratio shear viscosity 77 over entropy 
density s. Based on the correspondence between Anti-de- 
Sitter (AdS) space and conformal field theory (CFT), it 
has been conjectured [Sj] that all relativistic quantum field 
theories at finite temperature and zero chemical poten- 
tial have r]/s > To date, no physical system violating 
this bound has been found. 

Neglecting effects from bulk viscosity and heat con- 
ductivity, the energy momentum tensor for relativistic 
hydrodynamics in the presence of shear viscosity is 



T^- ^ {e+p)i 



(1) 



In Eq. ([T]), e and p denote the energy density and pres- 
sure, respectively, and u'^ is the fluid 4-velocity which 
obeys g^^u^u'^ — \ when contracted with the metric g^^- 
The shear tensor H^'' is symmetric, traceless (Ji^ = 0), 
and orthogonal to the fluid velocity, u^Ii^^ = 0. Conser- 
vation of the energy momentum tensor and equation of 
state provide five equations for the evolution of the 10 
independent components of e,p, u'', H^''. The remaining 
five equations for the evolution of H^" are not unambigu- 
ously agreed on at present [1, 0, S, S ■ The results in 



this work will be based on using the set of equations 
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-H''''[5i?lnr-V„u"], (2) 

where da is the covariant derivative, used to construct the 
time-like and space-like derivatives D = m"c?q and Vu = 



The remaining definitions are A 



U-) = V^w'^-f Vu^- IA^^Vqu" and the vorticity 
fj.L' = Vi,Wp — VpWi/. Both p and temperature T are re- 
lated to e via the QCD equation of state, for which we 
take the semi- realistic result from Ref. If the re- 

laxation time Tn is not too small, Eq. ([2|) are the most 
general shear viscous hydrodynamic equations that are 
causal and guarantee that entropy can never locally de- 
crease jl2j . Formally, Eq. ^ correspond to the relativis- 
tic Navier-Stokes equations in the limit rn — » 0, but con- 
tain corrections of higher order in gradients for rn > 0. 

Unfortunately, the initial conditions for a hydrody- 
namic description of an ultra-relativistic heavy-ion col- 
lision at RHIC are poorly known, so one has to resort 
to model studies. In order to describe Au-|-Au collisions 
at RHIC energies, one typically assumes the energy den- 
sity along the longitudinal direction (the beam-line) to 
be "boost-invariant" to first approximation [T^j- With 
this assumption, one still has to specify the energy den- 
sity distribution in the plane orthogonal to the beam line 
(the transverse plane). At present, there exist two main 
classes of models for this distribution, which we will refer 
to as Glauber-type and Color-Glass-Condensate (CGC)- 
type models. In the following, only Glauber-type models 
will be used. 

The Glauber-type models build upon the Woods- 
Saxon density distribution for nuclei, /ca(x) ~ !/[! + 
exp((|x|— i?o)/x)]i where for a gold nucleus {A— 197) we 
use i?o = 6.4 fm/c, x = 0.54 fm/c. Integrating the Woods- 
Saxon distribution over the longitudinal direction (corre- 



2 



1 ' 1 

r|/s=1.0 



• 1=1.5 fm/c 

♦ 'C=2fm/c 
■ T=3fm/c 




k [GeV] 



FIG. 1: Correlation function /(fc,T, to = Ifm/c) as a func- 
tion of momentum k, measured for our hydrodynamics code 
on a 64'^ lattice with a lattice spacing of lGeV~^ (symbols), 
compared to the "analytic" result from the linearized Eq. Q 
(full lines). The good overall agreement indicates the code 
is solving Eq. © correctly in the linear regime (see [3 for 
details) . 

spending to a Lorentz-contracted gold nucleus), one ob- 
tains the nuclear thickness function Ta(xj^). Contracting 
two Ta's for the collision of two gold nuclei at a certain 
impact parameter b, one can define number densities in 
the transverse plane, such as the density of participating 
nuclei npart , b) and the density of binary collisions 
ncoii(x_L,b) (see e.g. Ref. 



14l|). As an initial condition 



for hydrodynamics, one then assumes the energy density 
distribution e in the transverse plane to be proportional 
to either ripart or ncoii or a mixture of the two. In what 
follows, we will concentrate on the case e jt-coU, since 
for ideal hydrodynamics this provides a rough description 
of the centrality (or impact parameter) dependence of the 
total number of particles ("the multiplicity") measured 
at RHIC li*]. Finally, for VH one also has to provide 
an initial condition for H'^'^ . We choose the "minimal- 
ist assumption" 11^'^ = 0. While one realistically expects 
n'^'^ to be nonzero initially, this assumption translates 
to reducing the effect of viscosity, which can serve as a 
baseline for future studies. 

Because of boost-invariance, it is useful to work in the 
coordinates r = ^/t'^ — z^ and 77 — arctanh(z/t) rather 
than t, z. In these coordinates, boost-invariance dictates 
u'^ = 0, so because of u^^u'' = 1, the only non-trivial fluid 
velocities can be chosen as , , which are assumed to 
vanish initially. 

Before discussing results from the numerics, one can 
get some intuition of viscous effects on experimental ob- 
servables by imagining the system to have a friction force 
proportional to velocity. In a heavy-ion collision, the ex- 
pansion (at least initially) is strongest along the beam 
axis, therefore one expects viscosity to counteract this 
expansion. In r, 77 coordinates this is achieved by a reduc- 
tion of the effective longitudinal pressure p — IIJJ through 
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FIG. 2: Total multiplicity dN/dy and mean momentum for 
■K+,n-,K+,K-,p and p from PHENIX [H for Au-hAu col- 
lisions at ^/s = 200 GeV, compared to our hydrodynamic 
model for various viscosity ratios 77/s. 



n;j > 0. Since initially H; < IV^ but H*; = 0, the dif- 
ference between equilibrium pressure p and effective lon- 
gitudinal pressure has to appear as excess pressure in 
the transverse plane. Therefore, viscosity should lead to 
higher transverse velocities ("radial flow") as compared 
to ideal hydrodynamics, which is indeed the case 
Similarly, one can get an intuition of viscosity on elliptic 
flow V2 (the main angular modulation of radial flow for 
non-central collisions): having a stronger reduction effect 
on higher velocities, viscosity tends to decrease velocity 
differences and hence elliptic flow. This agrees with the 
qualitative trend found by Teaney 17| . 

To solve Eq. ^ and treat the freeze-out 



(see below), 

we have used a two-dimensional generalization of the al- 
gorithm outlined in Ref. ^1^. Details of the calculation 
will be given elsewhere [18j|. We have checked that our 
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algorithm agrees with the results from Ref. [19| for cen- 
tral collisions, when dropping the extra terms in Eg. (|2|). 
Also, our code passes the fluctuation test from Ref. \W^ . 
shown in Fig. [T] We thus have some confidence that our 
numerical algorithm solves Eq. ^ correctly. 

When solving the set of equations ([2]), we set the ratio 
r]/s to be constant throughout the evolution of the sys- 
tem, since modeling any space-time dependence would 
necessarily introduce more unknown parameters. There- 
fore, results on 77/s quoted below should be considered 
as mean values over the entire system evolution. 

To make contact with experiment, the hydrodynamic 
variables are translated into particle spectra via the 
Cooper-Frye freeze-out mechanism (adapted to VH 



L^ooper- 



see also [17|). For simplicity, we use a single 



freeze-out temperature Tf but include the effect of res- 
onance decays with masses up to 2 GeV on the spectra 



2IL 1221 . The normalization of the initial energy den- 
sity and Tf are chosen such that the experimental data 
on total multiplicity and mean transverse momentum 
< Pt > as a function of total number of participants 
.^Part = / c'^xj^ripart(x_L, b) are reasonably reproduced 
by our model (see Fig. [2]). We choose to fit to kaons 
rather than pions because the former are infiuenced less 
by Bose enhancement effects, which we have ignored [foj . 
Note that for simplicity our model does not include a 
finite baryon chemical potential, prohibiting us to dis- 
tinguish particles from anti-particles. As a consequence, 
results for protons cannot be expected to match exper- 
imental data. Starting from ideal hydrodynamics with 
a freeze-out temperature Tf = 150 MeV, we have found 
that reasonable fits to dN/dy and < pr > for VH can 
be accomplished by keeping Tf fixed and reducing the 
initial entropy density by 75 r]/s percent to correct for 
the viscous entropy production 

In Fig. [3] we compare our hydrodynamic model with 
the above fit parameters to experimental data on the in- 
tegrated and minimum bias elliptic flow V2, respectively. 
Shown are results for ideal hydrodynamics and VH for 
the initial condition e ~ ncoU at an initial time tq — 1 
fm/c. The results hardly change when assuming instead 
s ~ npart as initial condition (see also T^) or varying tq 
by a factor of two. Interestingly, we also flnd that chang- 
ing Tn hardly affects the results shown. Note that this 
depends on the presence of the terms in the last line of 
Eq. ([2]): if these terms are dropped, increasing rn tends 
to further suppress V2 in line with the trend found in [l^ . 

For the above initial conditions, we have noted that 
there is also hardly any effect from the vorticity term. 
This can be understood as follows: noting that for = 
the only non-trivial vorticity is w^*', which vanishes ini- 
tially because oi = — and forming the com- 
bination W^Du^ — W^Du^ we find -up to third order 



0.08 



0.06 



0.04 - 



0.02 




400 



Part 



25 



20 



c 15 



ideal 
Ti/s=0.03 
ri/s=0.08 
ri/s=0.16 
STAR 



ft 



10 




P^[GeV] 



FIG. 3: PHOBOS [2J] data on pr integrated V2 and STAR 
[2^ data on minimum bias V2, for charged particles in Au-f-Au 
collisions at -^s = 200 GeV, compared to our hydrodynamic 
model for various viscosity ratios 77/s. Error bars for PHO- 
BOS data show 90% confidence level systematic errors while 
for STAR only statistical errors are shown. 



corrections- 



Dp 

e + p 



= C'(H3). (3) 



This is the relativistic generalization of the vorticity 
equation, well known in atmospheric sciences [2^ . Start- 
ing from Lo'-^y ~ 0, Eq. ([3]) implies a very slow buildup of 
vorticity, explaining the tiny overall effect of the vorticity 
term in Eq. ^ . Note that upon dropping the assumption 



— 0, this term can become important 27[. 



From Fig. [3] it can be seen that the effect from viscos- 
ity on the elliptic fiow is strong, in line with estimates 
from Ref. [l3]. Data on integrated V2 is fairly well re- 
produced by a viscosity of "q/ s ^ 0.08 and - within sys- 
tematic errors - seems to be consistent with 77/5 ~ 0.16. 



4 



These values agree with recent estimates by other groups 

EE 



31|. How- 



|29|, l30| and a lattice QCD calculation 

ever, the comparison to data for minimum bias V2 in 
Fig. [3] suggests that the ratio of 77/s is actually smaller 
than the conjectured minimal bound 77/s = ^ ~ 0.08 As 
mentioned, this seems to be independent from whether 
one adopts rn = 6 77/(6 + p), the weak-coupling QCD 
result, or extrapolates to rn — > 0, which is very close 
to the AdS/CFT value found in Indeed, the mini- 
mum bias V2 seems to favor rj/s ~ 0.03, at least at low 
momenta, where hydrodynamics is supposed to be most 
applicable. Note that this result could change drastically 
if the minimum bias data were decreased by 20%, which 
is the estimated systematic error quoted in j25| . 

There are, however, a number of caveats that should 
be considered before taking the above numbers literally. 
Firstly, we have only considered Glauber-type initial con- 
ditions, and assumed 11^'' (tq) = 0. It has been suggested 
that CGC-type initial conditions lead to larger overall V2 
[33! which in turn would raise the allowed values for 77/s 
in our calculation. This is due to the larger eccentricities 
in this model [13] (note the issues raised in fsBl)- How- 
ever, larger eccentricities in general also lead to a faster 
build-up of transverse flow, which is further enhanced 
by viscosity. Thus, when required to fit all the data in 
Figs. [2] and [3l it is unclear whether this CGC-type model 
will predict substantially higher 77/s than found here. 

Secondly, we used VH until the last scattering instead 
of more sophisticated hydro-|-cascade models (e.g. [sl, 
[37II). We do expect changes in the extracted values of 77/s 
once a VH-|-cascade model description becomes available. 
Finally, at present we cannot exclude that effects not 
captured by hydrodynamics, such as strong mean-fields, 
distort our results. Work on QCD plasma instabilities 
and CGC dynamics might shed some light on this issue. 

To summarize, we have presented the first viscous hy- 
drodynamic fits to experimental data on the centrality 
dependence of dN/dy, < px > and V2 at top RHIC ener- 
gies. For Glauber-type initial conditions, we found that 
data seems to favor values for tj/s which are very small, 
below the AdS/CFT bound 0^. While suggested to be 
possible in [H, \^ , it will be interesting to see whether 
the above caveats - once addressed - can change our re- 
sults enough to accommodate viscosity equal or larger 
than the bound. In any case, we hope that our work can 
serve as a guideline to understanding the properties of 
the fluid created at RHIC. 
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